arXiv:1505.05634v4 [hep-ph] 31 May 2016 


RMS and charge radii in a potential model 


Tapashi Das 1 ’*, D K Choudhury 1,2 and K K Pathak 3 
department of Physics, Gauhati University, Guwahati-781 014, India. 

2 Physics Academy of North-East, Guwahati-781014, India. 

3 Department of Physics, Arya Vidyapeeth College, Guwahati-781016, India 

*email:t4tapashi@gmail.com 


Abstract 

The Dalgarno’s method of perturbation is used to solve the Schrodinger’s equation 
with the Cornell potential V ( r ) = — + br + c. The short range and long range effect of 

the potential is incorporated in the same wave function by using two scales r s and r L as 
an integration limit. With the wave function then the results for bounds on r.m.s. radii 
of various heavy flavored mesons are reported. We have also showed the relation between 
r.m.s. and charge radius of mesons. 
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1 Introduction 

Quantum Chromodynamics (QCD) [1] is established to be the successful theory of strong 
interaction and the remaining problem in onr understanding of QCD has to do with quark 
confinement in hadrons. It is expected that confinement is the dominant quark-antiqnark or 
quark-qnark interaction at large separation distances and therefore plays the vital role to explain 
the properties of highly excited (large sized) hadrons [2j. In the last few years, the heavy-heavy 
mesons (r/ c , T/A Xc, h c , i]b, Xb{ 1, 2, 3 P), hj,, 7 etc.) have been widely studied theoretically as well 
as experimentally. Again highly excited light-quark mesons show a hydrogen-like spectral pat¬ 
tern [3] and can be reproduced within a non relativistic constituent quark model framework 
which becomes asymptotically coulombic [!]. Thus the distance or scale behaviour of the QCD 
potential in studying the static and dynamic properties of hadrons is still a topic to explore. 
The method of perturbation in this context has the advantage of choosing parent (more domi¬ 
nant) and perturbative (less dominant) term from the QCD potential. The advantage of taking 
Cornell Potential for study is that it leads naturally to two choices of parent Hamiltonian, one 
based on the Coulomb part and the other on the linear term, which can be usefully compared 
[5J. Aitchison and Dudek in ref. [ 6 ] showed that with Coulombic part as parent (VIPT), bot- 
tomonium spectra are well explained than charmonium, where as charmonium states are well 
explained with linear part as parent. 

In this spirit a QCD potential model has been pursued for the two choices in seperation in 
ref. (7JEUH] by considering the Dalgarno’s method (TDi.Ilj of perturbation. The results include 
static and dynamical properties of heavy flavored mesons, like their form factors, masses, decay 
constants as well as Isgur wise function. Here in this work, we emphasise in the ground state of 
heavy-light mesons as well as cc and bb states, though cc and bb has numbers of excited states 


1 


for different T and ‘s’ values. 


However, while using perturbative method, one should consider two aspects of Quantum 
Mechanics: (a) The scale factor ‘c’ in the potential should not affect the wave function of the 
meson even while using perturbation theory to be compatible with quantum mechanical idea 
and 

(b) The specific choice of perturbative piece (coulomb or linear) should determine the pertur- 
batively compatible effective radial separation between the quark and the anti-quark. 

In the present work, we therefore show that only in the short distance range (0<r<r 5 ), linear 
potential is perturbatively compatible while for the alternative choice ( coulomb as perturba¬ 
tion), the corresponding range belongs to large distance {r L <r< oo). The exact magnitudes of 
r L and r s have explicit dependence on strong coupling constant a s , the confinement parameter 
b and on the scale factor c. 

The aim of the present paper is to outline the new features of the improved version of the 
model and present the prediction for the bounds on r.m.s. radii of pseudoscalar heavy flavored 
mesons (both heavy light and heavy heavy). Comparison is then made with those of other 
models available in literature. 

In section 2, we outline the formalism, while in section 3 summarize the results and discus¬ 
sions. Section 4 contains conclusions. 


2 Formalism 

2.1 Definition of r.m.s. radius 

The r.m.s. radius mm of the bound state of quark and anti quark like meson is defined as 



having radial wave function -0(r). 

2.2 Definition of charge radius 

The standard definition of charge radius [16] is 

(4) = —6^eF(Q 2 ) | 0 =,„ (2) 

where F(Q 2 ) is a meson form factor [TTl 118] . 

Using above equation, the charge radii of various heavy flavored mesons are calculated in 
ref. [15], m ■ 

2.3 Relation Between RMS and charge radius 

The charge radii can be measured with electromagnetic probe whether the RMS radii defined 
as in eqn.(l) cannot be measured with electromagnetic probe. The RMS radius is nearly the 
average <r 2 > of the quark wave function, which presumably may be determined in Quark Gluon 
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Plasma (QGP) experiments presently studied at LHC 11 9j . However a simple approximate 
relationship between the two can be found from the following equation 



<Vi> + 


4m? 


d A p | S(p) | 2 ( 


m. 


2/ 


derived by Godfrey and Isgur [2D]. 


(3) 


Here e* is the charge of the i th quark/antiquark, rq and rq are the distances of the two 
quarks/antiquarks measured from the centre of mass, m is the mass of the quark and E = 
a Jp 2 + m 2 . $(p) is the quark momentum distribution. The exponent / can be determined in a 
semi-empirical way. 


We obtain the inequality 

r 2 E > 6i<r 2 > (4) 

i 

making the transformation of the co-ordinate we choose one of the quarks/antiquarks located 
at the origin. It results in 

r|>e<r 2 > (5) 

where <r 2 > can be interpreted as the standard RMS as defined in eqn.(l) and e = 

It yields the desired results: the ratio of the charge radii and RMS radii of the meson should 
always be more than the square root of the charge of the quark/anti quark of the meson having 
the larger value. As an illustration for D + meson, the ratio should not be less than a/|. 

2.4 The total wave functions with an improved perturbative ap¬ 
proach 

The Schrodinger equation describing the quark-anti quark bound state is 

n 2 

-V 2 V>(r) + (E - VU(r) = 0 (6) 

2 m 

The standard QCD potential is defined as |16j 

Ary 

V (r) = ——^ + br + c (7) 

dr 

Where -| is due to the color factor, a s is the strong coupling constant, r is the inter quark 
distance, b is the confinement parameter (phenomenologically, b = 0.183GeR 2 |21j ) and £ c’ is 
a constant scale factor. 


For potential (7), we can make four choices: 

Choice-I: —as parent and br+c as perturbation 
Choice-II:br as parent and — 4^ + c as perturbation 
Choice-Ill: —+ c as parent and br as perturbation and 
Choice-IV: br+c as parent and as perturbation 
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It is well known in quantum mechanics that a constant term ‘c’ in the potential can at 
best shift the energy scale, but should not perturb the wave function. This important point 
was overlooked in earlier publications mm® of the subject, but the present work takes into 
account this. 


We use the Dalgarno’s method of perturbation to construct the wave functions. 
The wave function for choice-I [|7] including relativistic effect is 




Ah 


yfia, 


1 + cAq\/ Too - - pba 0 r 2 


e “o ; 


( 8 ) 


where A 0 is the undetermined co-efficient appearing from the series solution of Schrodinger’s 
equation and 


a o = 1 


-l 


h = 


777 . 1777.2 


(9) 


( 10 ) 


7771 + 7772 

Here t? 7 i and m 2 are the masses of quark and anti quark respectively and p is the reduced 
mass of the mesons and 


e=l-Wl- 


-cx< 


( 11 ) 


provides the relativistic effect 
ma. li za.tion constant given by 

Ah =- 


due to Dirac modihcation factor and TV] is the nor- 

1 ( 12 ) 


f 0 °°^r 1 + cA 0 V^o _ \^ba 0 r 2 


1 2 / \ —2e 2r 

— ) e a o dr 
ao J 


by 


Similarly, the wave function for choice-II upto 0(r 4 ) is calculated in ref. [23] and is given 


to^ (r) = — [l + Ai(r)r + A 2 (r)r 2 + A 3 (r)r 3 + H 4 (r)r 4 l + p 0 ] (—\ (13) 

r \a 0 J 

where A t [r] is the Airy function [2SJ and TV 2 is the normalization constant 

1 V 2 =---r- (14) 

/ 0 ro 4vr [1 + Ai(r)r + A 2 (r)r 2 + A 3 (r)r 3 + A 4 (r)r 4 ] 2 (Aifar + p 0 ]) 2 dr 

The co-efficients of the series solution 
the functions of a s , /i, b and c: 

Ai 

A _ -MW 1 - c) 

2 + Apik\ + p 3 k 2 


as occured in Dalgarno’s method of perturbation, are 


An = 0 


-2p^ 


2pih + pfk 2 


(15) 

(16) 
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A -2 / uW°A 1 
3 6 + 6p 1 k 1 + pfk 2 

—2pW°A 2 + 2 pbAi 
12 + Spiki + p\k 2 

The different parameters are given by, 

Pi = (2/x6)s 


(18) 

(19) 


( 20 ) 


Po 

(In our case n=l for ground state) 


37 r( 4 n — 1) 3 
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kl = i + k - 
r 


( 21 ) 


( 22 ) 


k 


0.3550281 - (0.2588194)p 0 
(0.2588194)p! 


ko = 


k 2 


(23) 

(24) 


W 1 

W° 



(25) 

(26) 


The above analysis illustrate that the scale factor ‘c’ has observable effect in the final wave 
function if we apply Dalgarno’s method of perturbation. Now in the present work we want 
to see whether the term ‘c’ has also its affect even if we consider ‘c’ in the parent part of the 
Hamiltonian. 


Using Dalgarno’s method of perturbation the total wave function for choice-III upto 0(r 4 ) is 


$ 


total . 


Ab 


r = 


yfna\ 


0 L 


1 - pib(i () r 2 - ^p 2 bca 0 r 4 


e a o 


where the normalization constant IV 3 is 
N 3 = - 


fo°° %[ l ~ \^a Q r 2 - ±p 2 bca 0 A] 2 e Z dr 

Similiarly, the wave function for choice-IV including relativistic effect considering upto 0(r 4 ) 


(27) 


(28) 


is 


i,btotai (r) = — [l + Ai(r)r + A 2 (r)r 2 + A 3 (r)r 3 + ff 4 (r)r 4 ] A^ppr + p 0 \ 
where iV 4 is the normalization constant given by 



(29) 
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iV 4 = 


f r long 4vr [1 + A\ 1 

(r)r + A 2 (r)r 2 + A 3 (r)r 3 + A 4 (r)r 4 ] 2 (A[pir + p 0 ]) 

dr 

The co-efficients are 


(30) 


O 

o 

(31) 


4 _ - 2 ^ 

2piAi + p\k 2 

(32) 


-2 pW 1 

2 + 4piAq + Pik 2 

(33) 


1 _ -2 p(W° - c)A 1 

3 Q + Qp x k 1 +p\k 2 

(34) 


—2 p(W° — c)A 2 + 2pbAi 

12 + Spiki + Pik 2 

(35) 


The equations (27) and (29) showed that even if we consider the scale factor ‘c’ in the parent 
Hamiltonian, then also it affects the total wave function. 

Therefore, Dalgarno’s method in general conflicts with the quantum mechanical notion that 
the scale factor ‘c’ should not have an observable effect except the energy shift. 

There are two possible ways of making the perturbative approach adequate to quantum 
mechanics. At phenomenological level one keeps the scale factor in the potential so that it can 
contribute to the short and long range of the inter quark distance without changing the total 
wave function in a true sense. A more theoretically satisfactory approach will be to limit the 
domain of the applicability of Dalgarno’s method, that it is not applicable to potential having 
non-zero scale factor. This approach is then corrected quantum mechanically. We take this 
latter view point in the present paper. 

For the validation of the quantum mechanical expectation we therefore consider c = 0 in 
the potential (7). In such situation choice-I & II and III & IV are identical. 

Now to find the cut offs we use perturbation conditions: 

Case-I: For coulomb as parent and linear as perturbation: 

- —>br (36) 

3 r 

and 


Case-II: For linear as parent and coulomb as perturbation: 


br> — 


4a* 
3r 


(37) 
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From (36) and (37) we can find the bounds on r npto which case-I and II are valid. Case-I 

S ^ vT an d case -H gives the cut off on the long 


gives the cut off on the short distance r^ ax < 
distance r£v > \ 


In table-1, we show the bounds on r s and r L in fermi which yields exact/most restrictive 
upper bounds of the quantities to be calculated. 


Table 1: r s and r L in fermi with b=0.183 GeV 2 


Scale 

r s = r L (fermi) 

c-scale 
(a s = 0.39) 

0.33207 

b-scale 
(a, = 0.22) 

0.24941 


Again the application of Airy function as meson wave function needs suitable cut off to 
make the analysis normalizable and convergent. We therefore set the cut off (r 0 ) in the range 
1 fin (5.076 Geld -1 ) |26j for our calculations. 

The equation (1) of section 2.1 is thus modified to 



3 Results and discussion 

Using equation (38) and with the help of the wave functions from equations (27) and (29) we 
compute the bounds on RMS radii of various mesons in the ground state (l = 0) and results 
are tabulated in table-2. However it is to be noted that r s and r L are the perturbative satura¬ 
tion lengths for coulomb parent and linear parent respectively. The proper perturbative range 
should be far less than the cut offs r s and r L . The upper bound corresponds to the sum of the 
contribution of the two options at short (coulomb parent) and long distances (linear parent), 
while the lower bound correspond to the minimum of the two (linear parent). 

The input parameters in the numerical calculations used are same as our previous work 
[271 [281 [29] : m u = 0.336GeU, = 0.336GeU, m s = 0.483GeU, m c = 1.55GeU and m;, = 

4.95GeU, b=0.183 GeV 2 . 
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Tabic 2: Bounds on RMS radii (■ r rms ) in fm with c=0 


(X. s 

Meson 

Bounds on 
{r r ms) in fm 

0.39 

7 r + (ud) 
K(us/ds ) 
D(cu/cd) 
D+(cs) 
3/4>(cc) 

0.2914-1.1334 

0.2940-1.1184 

0.3008-1.0924 

0.3091-1.0747 

0.3489-1.0530 

0.22 

Bt( be) 
r{bb) 

0.3841-1.0942 

0.4439-1.1331 


In tablc-3, we give the different model predictions of rms radii for heavy flavored mesons 
available in literature [301(311132]. 


Table 3: r rms values in fermi 


Meson 

r rms(f m ) 

Ref. [30] 

Ref. [31] 

Ref. [32] 

cc 

0.4490 

0.4453 

0.4839 

bb 

0.2249 

0.2211 

0.2672 


From tablc-3, we see that the value of rms radii for cc is well within the range of table-2, 
but the result for bb is out of range. A higher value of cut off parameter (ro) in this sense can 
improve the results for bb and then the other theoretical predictions can be accomodated within 
the range. 

In this work we have taken the cut off ro =lfm (5.076 GeV~ l ) [2_6J. However from the re¬ 
sult we can conclude that the scale for light-light, heavy-light and heavy-heavy mesons should 
not be same, rather it should be increased with the increase of the mass scale in a particular way. 

From table-2, the bounds on r rms for D + (cd) is 0.3008-1.0924 and hence the corresponding 
value of charge radii (te) should be more than this bound as total charge of the meson is +1 
and this fact is found to be within the prediction. Thus the validitity of equation (5) is also 
tested. 


4 Conclusions 

We have applied the perturbatively compatible approach to the QCD potential model studied 
earlier [3 El El EH EZ! and computed the upper and lower bounds of the RMS radii of heavy 
light mesons. The theoretical prediction of bound for cc meson is shown in table-3 . However 
the ‘c’ itself cannot be determined from the perturbative analysis. There is a suggestion in ref. 
[7] that the scale should not exceed ~ lGeV for heavy flavored mesons which is consistent with 
models of ref. [331 EH mi EH [32]. 

Let us now comment on the physical significance of RMS radii of heavy flavored mesons 
and their plausible experimentally observable possibilities. The present work will hopefully 




















be useful to extract information on RMS radii from the charge radii itself for heavy flavored 
mesons, where model predictions are available in current literature [351 ESI HH]- We also check 
the sensitivity of the cut off parameter (r 0 ) which is taken to be 1 fm in this work and the results 
show that the cut off parameter should not be same for light-light, heavy-light and heavy-heavy 
mesons, rather the cut off parameter should be increased from light to heavy mesons. Moreover 
it is to be noted that the coupling constant (ct s ) is chosen only in two scales: charmonium 
(a s = 0.39) and bottomonium (a s = 0.22) scales. A scale dependent a s [37] might be useful in 
this context. 

Although several theoretical models are available in literature, but there is no present exper¬ 
imental information either on charge radii or RMS radii of heavy flavored mesons. Hopefully, 
future experiments eg., PANDA [38], LHCb [19] will yield more information to test the validity 
of the simple relations as well as data. 
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